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Introduction. The development of a mathematical model for 
the groundwater level of a deposit of cement raw materials 
located in the Zadonian-Yelets aquifer, which is the principal 
domestic water supply source for the city of Lipetsk, is 
considered. Therefore, it is necessary to provide ongoing 
monitoring and to have the possibility to predict the water 
level under the field development. The work objectives are the 
identification and study of a dynamic neighborhood model 
with variable hierarchical neighborhoods of the groundwater 
level that enables to adequately predict value of the water level 
in the examined wells. 

Materials and Methods. The definition of a dynamic 
neighborhood model with variable hierarchical neighborhoods 
is given, differing by time-varying double-level neighborhood 
communications between the first- and second-level nodes. At 
each next discrete instant of time, the neighborhood model 
nodes change their state under the influence of the online 
parameters and node states included in their neighborhood. As 
a subcase, we consider a model with line state recalculation 
functions. Parametric identification of the dynamic 
neighborhood model consists in finding the system parameters 
for each second-level node, and is based on the ordinary least 
squares. 

Research Results. A linear dynamic neighborhood model with 
variable hierarchical neighborhoods for predicting the 
groundwater level in a cement raw material deposit located in 
the Zadonian-Yelets aquifer is developed. The software using 
C++ is developed for the parametric identification and 
simulation of the functioning of the dynamic neighborhood 
model under consideration. It enables to determine parameters 


of the node state recalculation functions for a given structure, 
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Beedenue. CTatba NocBaljeHa paspadoTKe MaTeMaTHYecKoli 
MOJeuIM YPOBHA MO J3CMHBIX BO, MeCTOPOXKeHHA LeEMeHTHOTO 
CbIPbA, PaCMONO%*KeHHOTO B 3aOHCKO-eJIeEHKOM BOJOHOCHOM 
TOpH30HTe, ABJIATOLHMCA OCHOBHBIM HCTOUHHKOM 
XO3AHCTBCHHO-IMTbeBOrO BOWOCHaOxeHHA Topoya JIumenka. 
Ilostomy Ha cTaqMH pa3pa0oTKH MecTOopoxyeHHa 
HeOOXOJHMO MpOBOZHTb NOCTOAHHbIM MOHHTOPHHT HM MMeTb 
BO3MO2%KHOCTb IIpOrHO3HpOBaHHA YPOBHA MO3eMHBIX BOL. 
Lemp padoTrr — upentTudpukalua u uccleqoBanne 
WMHaMHYeCKOH OKpecTHOCTHOM MoyemM Cc TepeMeHHbIMH 
WepapxXH¥ecKHMM OKPeCTHOCTAMH ypOBHA TMO3CMHBIX BOT, 
TO3BOJMIOWeH C BOCTATOUHOM TOUYHOCTbIO MpOrHo3HpOBaTb 
3HayeHHe ypOBHA BO B OOCJIeEAYeMbIX CKBaxKHHAX. 
Mamepuanst u  memoovi. IlpuseqeHo ompeyemeHue 
TMHaMHYeCcKOH OKpeCcTHOCTHOM MoyemM Cc TepeMeHHbIMH 
OTIMMYarIOMleHca 
BYXYPOBHeBLIMH 


MepapxuyecKHMu 
U3MCHAFOWUIMMUCA BO 


OKpecTHOCTAMH, 
BpemMeHH 
OKPeCTHOCTHBIMH CBA3AMH Mey y3laMH MepBoro HM BTOporo 
ypopHa. B kaxkybili clleqyroujMii MCKpeTHbIt MOMeHT 
BpeMeHH y3JIbl OKPeCTHOCTHOM MOeIH MeHAIOT cBOe 
cocTosHHe Mo BO3elCTBHeM TeKYIIMX ympaBieHHuii u 
COCTOAHHM y3I0B, BXOJAMJHX B HX OKpecTHocTH. B kayecTBe 
yacTHOrO cyly4ad paccCMOTpeHa MOjJeIb Cc JIMHeMHEIMH 
dyHKuMAmMH = =epecueta coctosHuli. [lapamerpwueckas 
WeHTHpUuKaLHA WHHAaMHMYeCKOM OKpecTHOCTHOH Moe 
3akIO4aeTCA B HaXOKZCHHH WapaMeTpoB CHCTeMBI IA 
K@kKQOrO y3la BTOporo ypoBHA HM OCHOBaHa Ha MeTOye 
HaMMeHbIUMX KBa[paTos. 
Pesyismamet —succiledoeanua.  Pa3padotaHa _—sIMHeliiHaa 
TMHaMH4ecKkad OKPeCTHOCTHad MOjeyIb Cc TepeMeHHbIMH 
WepapxXHuecKHMH OKPeCTHOCTAMH Id MporHosHpoBaHHua 
YPOBHA TNOZ3eMHBIX BOX MeCTOPOXKACHHA WeMeCHTHOTO CBIpbA, 
pacrosox*KeHHOTO B 3a]OHCKO-eJIeL]KOM BOJ]OHOCHOM TOpH30HTe. 
J\ia mapametpwyueckolt ueHTHpuKaywH UM MOjeIMpoBaHHsa 
(PyHKWMOHMpOBaHHA 
OKpecTHOcTHOH 


TWMHaMvyeckon 
TIporpaMMHoe 


paccmatTprBaemoli 
MOjeuIH pa3pa6oTaHo 
oOeceyeHve Ha sa3bIkKe C++, MO3BOUIAIOMIee JIA 3aaHHOl 
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and also to predict the model behavior in the operation 
process. A hierarchical structure is given, and a parametric 
identification of the linear dynamic neighborhood model of the 
groundwater level is carried out. After the parametric 
identification on the teaching data selection, the mathematical 
model is checked on the test sample. 

Discussion and Conclusions. The obtained average ratio errors 
of the identification and forecast suggest the developed model 
validity and enable to recommend it for predicting the 


underground water level of a cement raw materials deposit. 


Keywords: groundwater level, deposit of cement raw 
materials, dynamic neighborhood model with variable 
hierarchical neighborhoods, parametric identification. 
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CTIpyKTypbl HaxOQHTb apaMeTpbI (yHKMH mepecueta 
COCTOAHHH y3JI0B, a TakoKe TIPOTHO3HPOBaTb MOBeeHHe Moe 
B Tpoyecce dbyHkyHoHupoBanusa. Ipupeyena uepapxwueckad 
CTpykTypa HM UpoBeyeHa MapaMeTpueckad MeHTHpUuKalHA 
JmMHeHHOM = MHaMWYueCKOM OKpeCTHOCTHOH MOjleuIH ypOBHA 
TIOWBeMHBIx Boy. T[locne BbINONHeHHA TNapamMeTpHueckolt 
weHTupukaywH =69Ha =e OGyyalolyeli = BLIOOpKe =——-JaHHbIX 
MaTeMaTH4ecKasd MOJelIb IpoBepeHa Ha KOHTPOJIbHOM BEIOOpKe. 
O6cyacdenue  u Tlomyyenupie cpeqHHe 
OTHOCHTeIBHEIe § ONIMOKH HyeHTHpuKaWMH MM MporHno3a 
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Introduction. Zadonian-Yelets horizon is the main productive complex of limestone suitable for 
cement and metallurgical industries [1]. Limestone is light and light gray, with a yellowish tinge, medium-hard, 
fractured, fine-grained. You can often find hard limestone represented by silicified varieties. Limestone is often porous 
and cavernous. Clay rocks suitable for cement production lie among the quaternary deposits covering the limestone. 

Overburden rocks are represented by a fertile layer of soil, loam, substandard clay, sand and a layer of crushed 
stone in the upper part of the destroyed limestone. 

Hydrogeological conditions of the deposit are simple. Sokolsco-Sitovsky deposit is in in geomorphological 
terms in the slope part of the valley of the Voronezh River, which determines the hydrogeological situation in the field 
area [2]. 

The aquifer is fed by infiltration of atmospheric precipitation and by absorption of flood flow through the 
valleys of gullies and ravines. Groundwater discharge is in the Voronezh River. 

The current position of the groundwater level was studied by monitoring observations in 7 wells evenly located 
on the field area. The level of groundwater in the low-water period varies from 102.5 m to 109.7 m. 

Currently, the Zadonian-Yelets aquifer is the main source of the domestic and drinking water supply of the 
Lipetsk city. In this regard, the development of limestone can be made only in the water-free part of them with high 
pillar of at least 2.0 m, which is recommended by the sanitary service of the Lipetsk region. 

Water inflow into the developed pit of the deposit occurs only due to atmospheric precipitation. In this regard, 
at the stage of development of the field, it is necessary to constantly monitor the state of groundwater of the Zadonian- 
Yelets aquifer and to be able to predict the groundwater level. 

Dynamic neighborhood models with variable hierarchical neighborhoods are used in this paper to model the 
groundwater level of the field under consideration. These models allow modeling complex spatially distributed 
processes and objects that change their state in time [3-8]. 

Neighborhood models were first proposed in the late 90-ies of XX century [9]. The basic definitions and 
algorithms of the theory of neighborhood modeling are given in [9-12]. The concepts of “neighborhood”, 
“neighborhood communications” are considered in [13-14]. They use agents that move around the neighborhood and 
interact with each other according to certain rules. Today, the theory of neighborhood modeling is being actively 
developed. There appeared dynamic [6-8], non-deterministic models [11], models with variable neighborhoods [15]. 

Materials and Methods. In this paper, we use dynamic neighborhood models to predict the state of spatially 
distributed systems. In them, each node is an independent object, functioning in time and related by neighborhood 
communications to the other objects of the system. This distinguishes the models under consideration from the widely 
used today neural networks, which can be used to simulate the operation of each object or node separately. The method 
of neighborhood modeling is intended for simultaneous joint modeling and prediction of the behavior of all elements of 
a distributed system. 
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Next, we consider the dynamic neighborhood models “input-state” with variable hierarchical neighborhoods. 
They differ in time-varying double-level neighborhood communications between the first- and second-level nodes, and 
in line state recalculation functions. They enable, in comparison with the known one-level neighborhood models, to 
perform forecasting with higher accuracy. 

The dynamic neighborhood “input-state” model with variable [15] hierarchical neighborhoods can be specified 
by the set NS, =(N,X,V,G, X[0],t) , where: 

1) N =(4,0,,0,,0,,) is a two-level structure of the neighborhood model; A= {a,,a,,...,4,} is a set of nodes 
of the first level; O, and O, are neighborhoods of communications of nodes by states and by controls, respectively; 


O,,,. are hierarchical neighborhood communications between nodes. Each node a, ¢A has its own neighborhood 


defined by states O,[a,] CA and controls O [a,]c A; O, = UO, (a), O.= YO, [a,]. 


Some nodes of the first level a, € A are assigned a set of nodes of the second level O,,,[a,] = {a;,...,a°}. All 


nodes a, €O,[a;], a, €O,[a,] have an impact on the second level nodes a €0.fe), 
At each point in time, only a single active node is defined b =1,...,c such that a? € O,,,[a,,¢]. 
The node a? is active at the time f¢ if it meets the specified activation condition Baa tal = true. For all nodes of 


the second level a? O,[a?]=O,[a,]; O,[a?]=O,[a, ]. 


Yr 
2) XeR" is a block vector of states of the neighborhood model in real time, each block of which 


X[a,|= X[i]e R” is a vector of states in the node a, of the system i =1,...,n. 


ym 
3) VeR™ isa block vector of controls in real time, each block of V[a,]=V[i]¢ R” is a vector of controls 


in the node a, ofthe system i=],...,.n. 


4) G:X, xV, +X isa vector function of states recalculating of the neighborhood model, where X, is a 
set of states of the nodes of the first level, included in the neighborhood O, ; V,_ is a set of controls of the nodes of the 


first level, included in the neighborhood O, . 


For the nodes of the first level a, € A, the function G, will be: 
X[t+1,=Gld =>) X[t+1,77]= > GA, (1) 
b=1 b=1 


where G?:X,, OF ee [i’] is a state recalculation function for the second-level node a?. 


For each node of the second level a? € O,,[a,,¢], the function G? in real time in the linear case has the form: 


teri 


Xe+h’]= >) etl? JAB /+ >) ell? AV A+ 220"), (2) 


a;€Oy [ay ] 4 €0, [ay] 
where a,,a, € A (ik = Length) are nodes of the first level of the model; X[t,i’]¢R” is the state in the node 


a’ at time ¢; V[t,i’]€ R” is the control in the node a? at time ¢; g°[i’,jJeR”””’, g?[i’,k]eR™™, gf’ ]eR’™ are 
model matrix-parameters. 
For each second-level node a? ¢ O.,,[a,,t], the function G? is currently zero, that is X[¢+1,7’]=G?[¢]=0. 


jeri 


va 
5) X[0]<¢R“ 1s the initial state of the model. 


6) ¢ is the current discrete time of the model. 

The structure of the neighborhood model can be represented as a single two-level graph of the structure of the 
neighborhood model-oriented graph with two types of arcs: state and control actions — or two-oriented graphs — 
external and internal structures. 

Parametric identification of the dynamic neighborhood model [16-17] is in finding the system parameters for 
each node of the second level, and it is based on the least squares method: 


E=>)|Xt+1i)-GJs]] > min. 
i=l 
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Research Results. Let us consider the linear dynamic neighborhood model of the level of cement raw 
materials underground waters deposits. Fig. 1 shows the graph of the external structure of the model. Node a, is 
external environment. Currently, 7 wells are used for groundwater level monitoring purposes, which in Fig. 1 
correspond to the nodes a,—a, . The control actions of the nodes V[t,7]<¢ R'° consist of the amount of precipitation in 


mm and of the average daily air temperature in °C for the last 5 days before the measurement of the water level in the 
wells at time ¢. The states of the nodes X[t,i] € R are equal to the water level in the well i at time? , 7 =2,...,8. 





Fig. 1. Graph of the outer structure of the neighborhood model 


[4,]={a',a?} and 


correspond to the positive and negative average daily air temperature. Their structure is shown in Fig. 2. 


V{i| 


The first level nodes of the neighborhood model a, (i=2,...,8) are hierarchical: O 


ier 





Fig. 2. Graph of the inner structure of the first level nodes 


The system (2) for each node of the second level a? of the surrounding groundwater level model will have the 
following form: 
XUe+L?]= gh P+ lh Wb?) + 2), (3) 
where g°[i,i]—eR; g°[i,i]eR”"; ge? [i]eR; i=2,...,8; b=1,2. 
The program in C++is developed for the parametric identification and performance simulation of the 
considered dynamic neighborhood model. The initial data for identification are the structure and the training sample. 


The program allows you to find parameters of the recalculation functions of node states, as well as to predict the 
behavior of the mathematical model in the process of operation. 
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After the parametric identification on the training data sample, the obtained model was tested on the control 
sample. The normalized initial and model values of the groundwater level for the node a, for training and control 


sampling are shown, respectively, in Fig. 3 and Fig. 4. 


=—@=X[2] —=X"(2 


Fig. 3. Normalized source and model values of the groundwater level for the node a, 


=—X[2] = x*[2 


Fig. 4. Normalized source and predicted values of the groundwater level for the node a, 


The average relative error of identification (prediction) of the neighborhood model is calculated by the 
formula: 
MX [t+Li]-X,[t+Li 
ae 1 >» mf + Li] mie 4! 100%, 
Mn X,[t+Li] 


k=1 i=l 





where X,,[f+1,i] is the state of the node in the m-th training (control) sample tuple; Xx [iri i] are model 
values of the node state a,; M is the volume of the training (control) sample. 


The average relative error of prediction is shown in Fig. 5. 
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Fig. 5. Average relative prediction error for each tuple of the control sampling 


The average relative error of identification was 0.19%, the average relative error of prediction was 0.23%, 
which indicates the adequacy of the developed model and allows us to recommend it for predicting the level of 
underground waters of the cement raw material field. 

Conclusion. A linear dynamic neighborhood model with variable hierarchical surroundings to predict 
groundwater level deposits of cement raw materials, located in Zadonian-Yelets aquifer, is developed. 

The software in C++programming language is developed for the parametric identification and performance 
simulation of the considered dynamic neighborhood model. 

The implemented model predicts the level of groundwater in the analyzed wells with sufficient accuracy and 
can be effectively used to predict the level of groundwater deposits of cement raw materials. 
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